Intrinsically s-wave like property of triplet superconductors with spin-orbit coupling 
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We have studied a general property of unconventional superconductors with spin-orbit coupling by 
solving the Bogoliubov-de Gennes equation and found that the two of four eigenfunctions for triplet 
pairing coincide with those for s-wave pairing when a relation holds between pairing symmetry and 
spin-orbit coupling, indicative of an intrinsically s-wave like property of triplet superconductors. 
Applying the result, we have studied the tunneling conductance in normal metal / insulator / 
unconventional superconductor junctions with a candidate of pairing symmetry of CePtsSi, d (k) oc 
yiky—yk^. The effect of Rashba spin-orbit coupling (RSOC) is taken into account. We have found 
that the RSOC induces a peak at the energy gap like s-wave superconductors in the tunneling 
spectrum, which stems from the intrinsically s-wave like property of triplet component of CePtsSi. 
As a result, the tunneling spectrum has s- and p-wave like features. This may serve as a tool for 
determing the pairing symmetry of CePtsSi. 



The recent discovery of CePtaSi has attracted much 
attention because it is the first heavy ferniion super- 
conductor without inversion symmetry^. It is predicted 
that this effect causes the Rashba type spin-orbit cou- 
pling (RS0C)2i^. Theoretically suprconductivity under 
broken mirror symmetry shows some novel phenomena 
due to the RS0C45.. Howewer the effect of the other 
types of spin-orbit coupUng (SOC), e.g., the Dressel- 
haus type SOC®, in superconductors without inversion 
symmetry is not studied well. It is important to study 
the general case because the pairing symmetry and the 
type of SOC depend on the crystal structure. Frigeri 
et al. showed that the pairing symmetry of CePtsSi is 
given by d (k) oc Stky—yk^^ where x and y are the unit 
vectors. This triplet pairing symmetry is, however, not 
conclusive^. It is desirable to calculate the physical quan- 
tities of CePtsSi with candidates of its pairing symmetry 
and compare them with experimental results. 

All these motivate us to study the general property of 
unconventional superconductors (USs) with SOC. After 
we clarify it, we apply the result to the tunneling junc- 
tions with CePtsSi and find a peak structure at the en- 
ergy gap in tunneling spectrum even for the triplet pair- 
ing. This may serve as a tool for determing the pairing 
symmetry of CePtsSi because the tunneling spectroscopy 
in superconducting junctions is a powerful method to 
study its pairing symmetry. 

In normal metal / supercunductor (N/S) junctions 
Andreev reflection (AR)* is one of the most important 
process for low energy transport. Taking into account 
the AR, Blonder, Tinkham and Klapwijk (BTK) pro- 
posed the formula for the calculation of the tunneling 
conductance^. This method makes it possible to clarify 
the energy gap profile of superconductors. It was ex- 
tended to normal metal / unconventional superconduc- 
tor (N/US) junctionsiSiii in order to study the properties 
of unconventional superconductors. In fact, calculation 
or measurement of the tunneling conductance in N/US 
junctions are useful to study the symmetry of the pair 
potential of USs because the tunneling conductance is 



sensitive to the pairing symmetry due to the formation 
of Andreev resonant statesiSiii. However, there is no 
theory considering the effect of the SOC in the tunnel- 
ing junctions. In this paper we generalize the theory 
in RefjiSiii by incorporating the RSOC and apply it to 
junctions with a new superconductor CePtsSi using the 
triplet pairing symmetry (d (k) cx xky—ykx). 

We will first clarify a general property of USs with SOC 
by solving the Bogoliubov-de Gennes (BdG) equation. 
Applying the result, we calculate the tunneling conduc- 
tance in normal metal / insulator / unconventional su- 
perconductor CePtsSi junctions in the presence of RSOC. 
We find that the RSOC induces a peak at the energy gap 
in the tunneling conductance even for the triplet pairing. 
The result can be explained by the general property of US 
with SOC. Futher we estimate parameter values for the 
tunneling junctions with CePtaSi. The present results 
may give useful information for the analysis in exper- 
iments to determine the pairing symmetry of CePtaSi. 
Although we use a simple model for CePtsSi which has 
a complicated band structureiS, we believe our results 
grasp the essence of the physics. Below we focus on the 
zero temperature regime. 

Let us start studying a general property of USs with 
SOC. Consider an effective Hamiltonian for BdG equa- 
tion with SOC. The Hamiltonian reads 



H 



H{k) A(k) 
-A*(-k) -H*{-k) 



(1) 



with H (k) = ^fc + V (k) -cr , and A (k) = iAay for singlet 
pairing or A (k) = (d (k) • a) iuy for triplet pairing. Here 
^fe, k and a denote electron band energy measured from 
the Fermi energy, electron momentum and Pauli matri- 
ces respectively. The second term of i?(k), V (k) • cr, 
represents the SOC. For example, V (k) cx xky—yk^ cor- 
responds to the RSOC and V (k) cx xk^—yky does the 
Dresselhaus type SOC. In this work we focus on the uni- 
tary states. We assume d (k) || V (k) where both vec- 
tors have only real number components and V (— k) = 
—V (k) which breaks inversion symmetry but conserves 



2 



time reversal symmetry. The condition d (k) || V (k) 
gives the highest Tc of the US'^ . The BdG equation reads 



H 



v± 



E4 



v± 



for clcctron-hke quasiparticles, and 

H ( "y^.^^y \ = -E± 

for hole-hke quasiparticles, 



(TyV±ay 

ayU±ay 



(2) 



(3) 



(a-±|V(k)|)^ + |Ar, u± = u 
and v± = (l±V(k) -cr 



with E± 

1 ± V (k) • 



1(1 + ^ 




Here 



1 - 



E± 



V(k) = V(k)/|V(k)| and |A|^ = iTrAAt. This 
shows that there exist the independent four eigenfunc- 
tions: electron- and hole- like quasiparticles with the 
eigenvalues E±. 

Let us discuss the property of the four eigenstates. For 
singlet pairing, we can find 



v± 



At / - 
^^o^(l±V(k).a 



'iayVQ ^1 ± V (k) • a 



with A = iAffy while for triplet pairing we get 



(4) 



v± = -iayvo [Y (k) ■aj[l±Y (k) • aj 

= T^^yvo (1 ± V (k) • a) (5) 

with A = (d (k) • a) iuy. Thus we can find that the eigen- 
functions with have the same form for singlet and 
triplet pairings. This indicates that the two of four eigen- 
functions for triplet pairing coincide with those for singlet 
pairing when the magnitude of the gap, |A|, has the same 
dependence on k, especially with those for s-wave pair- 
ing when |A| is independent of k. This feature gives an 
intrinsically s-wave like property of triplet superconduc- 
tors with SOC. As an example of its manifestation, we 
will calculate the tunneling conductance in the normal 
metal / unconventional superconductor (CePtaSi) junc- 
tions because CePtaSi is considered to satisfy the above 
conditions to have an s-wave like property^. 

We consider a two dimensional ballistic N/US junc- 
tions. The band structrure of CePtaSi shows that the 
main contribution to the density of states stems from 
the (3 band as shown in Refii^. The /3 band has a three 
dimensional complicated structure. However, its volume 
is large for large kz and hence the most important part 
is the (i band for large where the dependence of the 
(3 band on kz is weak. Therefore we focus on this part 
and assume the two dimensional N/US junctions as a 
first step. The N/US interface located at x = (along 



the y-axis) has an infinitely narrow insulating barrier 
described by the delta function U{x) — U5{x). We 
choose d (k) = 1^ {Siky—ykx) and V (k) = A (xfcy— yfc^;) 
with Rashba coupling constant A. The eigenf unctions of 
the Hamiltonian are ^ (uq, — iaj~^uo, iai^vo, vq), 

T ii^. „, „, „, „, \ T , 



with ai(2) = uo = 



(ia2Wo, -I'D, "0, «a2Uo) 



,+ 



Wo 



and A:i(2)± = fci(2)e**''i<2) . Here we put i?+ = 
^1(2) is an angle of the wave with wave number fci(2) 
with respect to the interface normal, fcj? is Fermi wave 
number, and m is effective mass in US. Velocity operator 
in the x-direction is defined asiSi — . 

Wave function ip{x) for a; < (N region) is represented 

as 




for an injection wave in up (down) spin states with 
kpx = kp cos 9x where 9 is an angle of the wave with 
wave number kp with respect to the interface normal in 
the N region, a and b are AR coefficients, c and d are 
normal reflection (NR) coefficients. Because we focus on 
the low energy transport compared to the Fermi energy, 
we neglect the difference of the wave number between 
electron and hole. 

Similarly for 2: > (US region) ip{x) is given by the 
linear conbination of the eigenfunctions. Note that since 
the translational symmetry holds for the y-direction, the 
momenta parallel to the interface are conserved: ky = 
kp sin 9 = ki sin 9i — k2 sin 02- 

The wave function follows the boundary conditionsi^: 



V'(a;)L=+o = V'(2^)L=_o 

10 
h 2mU ( 1 



0-10 
-1 



V'(o). 



(7) 



(8) 



Applying the BTK theory to our calculation, we obtain 
the dimensionless conductance represented in the form: 



E 



l + |a|V|6|'-|c|^-|d|^ 



(9) 



We define the normalized conductance as ctt = crs/ <Jn 
where ctat is given by the conductance for normal states. 
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and parameters as [3 



2m X 



and 



i.e., as for A 

Z = 9 is the effective mass in N divided by that in 

US. 

In the following we study the normalized tunneling 
conduntace ctt as a function of bias voltage V . For 
Z = 10 and g — 0.1, the magnitude of ax at eV — A 
increases as increasing /3 and finally a peak structure ap- 
pears (Fig. n(a)). On the other hand we can't find such 
a peak for Z = and g ^ 0.1 (Fig. [11(b)). Next we 
will study the case of a larger effective mass in US to 
check its effect on the conductance. Figure |21 shows the 
conductance for (a) Z = 10 and g = 0.01, and (b) Z = 
and g = 0.01. In both cases a peak emerges at eV = A 
and becomes sharpe for large f3. 
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FIG. 1: (color online) Normalized tunneling conductance with 
g = 0.1 for (a) Z = 10 and (b) Z = 0. 

Let us explain the origin of the peak at eV = A in 
Figs, n and |21 A coherent peak at eV = A in the 
tunneling conductanace also appears in N/S junctionsSr 
even in the presence of RSOC^''. The two of four 
eigenfunctions of the BdG equation with RSOC for the 
pair potential we consider coincide with those for s-wave 
pair potential. We can choose the electron-like and 
hole-like quasiparticle states with wave number ^2 as the 
two eigenfunctions. Acutually, using the expressions of 
u± and v± , we get the eigenfunctions of the Hamiltonian 
for s-wave pairing: {uq, — «aj~^wo, ia^^vo, vq), 

^(iaiuo, vo, -uo, iaiuo), ^(102^0, -"o, uq, ia2Uo). 
This implies that US with RSOC has an intrinsically 
s-wave like property. As /3 increases, the contribution 
of the eigenstates with wave number k2 is much more 
dominant than that with wave number ki because ^2 
(fci) is an increasing (a decreasing) function of /3. There- 
fore the behavior of the conductance becomes similar to 
the one for the s-wave pairing with increasing RSOC in 
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FIG. 2: (color online) Normalized tunneling conductance with 
g = 0.01 for (a) Z = 10 and (b) Z ^ 0. 



US. Note that near zero voltage the mig gap Andreev 
resonant states greatly change the spectrum compared 
to the one for the s-wave junctionsiSiii. Thus there 
is a qualitative difference in the conductance between 
two pairing states near zero voltage. For understanding 
the underlying physics, it is useful to check the case 
of very large RSOC, though it may be unphysical. 
Figure |3I displays the tunneling conductance with the 
same parameters in Fig. ^ (a) with very large RSOC. 
Apparently an s-wave-like tunneling spectrum emereges 
as increasing /3, which confirms the above explanation. 
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FIG. 3: (color online) Normalized tunneling conductance for 
Z = 10 and £/ = 0.1. 



Note that this feature is unique to the pair potential 
considered. As a reference, we will consider the case of 
other pairing symmetries with the same parameters as in 
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Fig. n] (a) . We choose 
in Fig. 0] (a) and 

in Fig. 21 (b). In both cases there is no qualitative 
change by the RSOC in the tunnehng conduntance. 
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FIG. 4: (color online) Normalized tunneling conductance for 
other pairing symmetries with Z = 10 and g = 0.1. (a) 
d (k) = -|| fc.y and (b) d (k) = -^fe,y. 

Now let us examine the parameter values suitable to 
CePtsSi. In the case of CePtsSi, we estimate g ~ 0.01 
from the specific-heat measurement and P ^ 0.3 from 



the band calculationia. The corresponding results are 
shown by solid curves in Fig. [51 (a) and (b), where the 
conductance clearly reflects the features of s- and p-wave 
supercondoctors: a broad peak around zero voltage and 
a peak at the energy gap. Recent NMR study of CePtsSi 
showed that 1/TiT of CePtsSi exhibits a small peak 
just below Tc~^^ like the Hebel-Slichter peak in s-wave 
superconductorsii. This similarity between CePtaSi and 
s-wave superconductors is consistent with our results. 

In summary we solved the BdG equation for US with 
SOC and found that the two of four eigenfunctions for 
triplet pairing have the same form as those for singlet 
pairing under the following conditions: (i) d (k) and 

V (k) have only real components, (ii) d (k) || V (k). (iii) 

V (— k) = —V (k). Moreover, if (iv) the magnitude of the 
gap, |A|, is independent of k, the two of four eigenfunc- 
tions for triplet pairing coincide with those for s-wave 
pairing. This feature gives an intrinsically s-wave like 
property of triplet superconductors with SOC. Note that 
this is essentially different from the recent work where the 
mixture of singlet and triplet pairing is discussed^. We 
found that triplet superconductor itself has a singlet-like 
property. As an example of its manifestation, we stud- 
ied the tunneling conductance in N/US junctions with 
RSOC and a candidate of pairing symmetry of CePtaSi. 
This compound is considered to satisfy the above four 
conditions^. We found that the RSOC induces a peak at 
the energy gap in the tunneling conductance as found in 
s-wave superconductor junctions, which results in the s- 
and p-wave like features in the tunneling spectrum. This 
may serve as a tool for determing the pairing symmetry 
of CePtaSi. 
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